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p I Abstract 

D 
^S^ We perform a next-to-next-to-leading order (NNLO) QCD analysis of the charni- 

top-quark contribution rjd to the effective \AS\ = 2 Hamiltonian in the Standard 

CN Model, rjct represents an important part of the short distance contribution to the 

^ parameter ex ■ We calculate the three- loop anomalous dimension of the leading op- 

QQ erator Qs2, the three- loop mixing of the current-current and penguin operators into 

^O Qs2, and the corresponding two-loop matching conditions at the electroweak, the 

^-% bottom-quark, and the charm-quark scale. As our hnal numerical result we obtain 

^^ Vet = 0.496 lb 0.047, which is roughly 7% larger than the next-to-leading-order (NLO) 

§ value r?^^^ = 0.457±0.073. This resuhs in a prediction for \eK\ = (1.90±0.26) x 10-^ 

y—i which corresponds to an enhancement of approximately 3% with respect to the value 

j> obtained using ?7^'^^. 

X 

c^ 1 Introduction 

Indirect CP violation in the neutral Kaon system was discovered by Christenson, Cronin, 
Fitch and Turlay in 1964, who observed the decay of a K^ into two pions fl]. This decay 
would be forbidden in the case of exact CP symmetry. The parameter ex measures indirect 
CP violation and is defined by 

^K = ^r7T? -7— T V (1-1) 
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Figure 1: The AS" = 2 box-type diagram with internal up, charm, and top contributions is 
expressed as a sum of box-type diagrams proportional to Af , A^, and At Ac, respectively, using the 
GIM mechanism. 



via the ratio of the respective decay amphtudes of a K^ and a Ks decaying into a two- 
pion state of isospin zero in such a way that direct CP violation is absent to a good 
approximation. 

The parameter ex is measured with high accuracy: The value quoted by the Particle 

(43.5±o.7)° ^_ Whereas until about a decade 



Data Group is ex = (2.228 ± 0.011) x lO'^ x e' ^ 

ago the numerical value of ex was used as an input to determine the Standard Model 
parameters, nowadays it plays a central role in constraining models of new physics: The 
near diagonality of the Cabibbo-Kobayashi-Maskawa (CKM) matrix leads to a suppression 
in the Standard Model, while ex can be predicted very reliably. 

For the theoretical prediction it is useful to express ex in terms of {K^\'Hr^~ \K^) = 
2MxMl2i the matrix element of the AS* = 2 effective Hamiltonian, and write: 



ex 



e^'^^sinc 



/Im(M^ 
V AM, ^ 



K 



:i.2) 



Here Mx is the neutral Kaon mass and AMx the Kaon mass difference, the phase of ex 
is 0e = 43.5(7)° [^ and ^ = ImAo/ReAo ~ is the imaginary part divided by the real 



part of the isospin zero amplitude Aq = A{Ks — )• (7r7r)/=o) 
45°,^ = 0)1 encompasses the change of \ex\ if the values 



The ratio k^ = 
k = 45° and ^ 



\,SM 



/ex{,<Pe = 

are used 



in (1.2), as has been done in most of the older analyses, instead of the exact values. The 
authors of Reference p5| give the value of k^ = 0.94 ±0.02 in the Standard Model, including 
in their analysis also contributions of higher-dimensional operators to the absorptive and 
dispersive part of the K^ - K^ mixing amplitude. 

The box diagram of Fig. [l^ gives the leading contribution to the effective Hamiltonian 



n 



|A5|=2 
/=3 



and the parameter M12. It is proportional to a sum of loop functions times CKM 



factors, which, using Xi = V*.Vid and Xi = mf/My^ we can write as: 

^ \u^\u^S{Xu,,Xu^)=:\lS{xt) + \lS{Xc)+2\c\tS{xt,Xc), (1.3) 

where S denotes the contributions of the individual box diagrams. After the Glashow- 
Ihopoulos-Maiani (GIM) mechanism has been used to ehminate A^ = — A^ — Ac it comprises 
the top-quark contribution - proportional to A^ (Fig. [ip), the charm-quark contribution - 
proportional to A^ (Fig. [ip), and the charm-top-quark contribution (Fig. [l|i) - proportional 
to AcAj. The resulting loop functions S'(xi, Xj) = S{xi, Xj) — S{xi, 0) — 5(0, Xj) + S'(0, 0) and 
S{xi) = S{xi, Xi) are suppressed by the smallness of the quark mass rrii if x, is significantly 
smaller than one. This, together with the severe Cabibbo suppression of the CP violating 
top-quark contribution, lets all three contributions compete in size for ex'- 

Im (A?5(xt) + XlS{x,) + 2XAcS{xt,x,)) 

(1.4) 



C(.VC.(.= 3;)h-o(.«^,o.(^)) 



where A = |\4s| ~ 0.2255. The diagram of Figure [l^ induces a large logarithm logrric/Mw 
only for the charm-top-quark contribution: the large logarithm from the up quarks in 
Fig. uh is power suppressed by Aq(^j-,/M^, while the GIM mechanism cancels a potential 
log rric/Mw between the diagrams with both one up and one charm quark and the diagram 
with only internal charm quarks. 

This can be reformulated in the language of an effective theory: the dimension-six 
penguin as well as the current-current operators, which have tree-level Wilson coefficients, 
mix only into the charm-top-quark contribution, via the bilocal mixing in Fig. [2^, yet do 
not induce large logarithms times tree- level Wilson coefficients proportional to A^ and A^. 
QCD corrections do not change this picture but only induce the well known renormalisation 
group effects for the AS" = 1 effective Hamiltonian kj and for the AS* = 2 Operator Qs2 
(Fig. |2b). A leading order (LO) analysis of the charm-quark and top-quark contribution to 
e/c then requires a one-loop calculation both for the matching at Hw and for the running, 
and for the charm-quark contribution also for the matching at f^c (Fig. [2k). This is in 
contrast to the charm-top-quark contribution where a tree-level matching at fiw and /Xc is 
sufficient at LO. 

After integrating out the charm quark the AS = 2 effective Hamiltonian reads 

^?=3=' = ^M^ [>^lVccS{x,) + \httS{xt) + 2K\t7],tS{x,, xt)] 6(/i)Q52 + H.c. + . . . (1.5) 
where Gp is the Fermi constant and 

QS2 = {-SLlM ® {'SLl^dL) (1.6) 

is the leading local four-quark operator that induces the \AS\ = 2 transition, defined in 
terms of the left-handed s- and (i-quark fields. The QCD and logarithmic corrections are 






a) b) 

Figure 2: Dimension-six current-current and penguin operators mix at LO into Qs2 with a CKM 
factor AiAc in a). Integrating out the charm quark results in similar diagrams for the LO and 
NLO matching of the contribution proportional to XtXc and A^, respectively. A sample diagram 
which is relevant to the LO evolution of Qs2 is shown in b). 

known at LO [s] and NLO and are parameterised by r]cc = 1.43(23) [oj, 7]ct = 0.47(4)[^ [t], 
and rftt = 0.5765(65) fs]. The parameter 6(/i) is factored out such that the bag factor 

is a renornialisation-group invariant quantity, which can be calculated on the lattice with 
high precision - see for instance [9]. Here fx is the Kaon decay constant. 

Finally note that 'Hf£^'^ also contains higher-dimensional operators and current-current 



operators with up quarks, as indicated by the ellipses in Equation (1.5). At LO in the 1/Nc 
expansion {N^ being the number of colours) only one higher-dimensional operator is present 
and its matrix element is estimated in |3l[l0] to result in a 0.5% enhancement of e^. 

In view of the improvements on the long distance corrections achieved in the recent 
years, the short distance contributions should be reconsidered. In this work we calculate 
the NNLO corrections to the charm-top contribution rjct- The NNLO corrections to the 
charm-quark contribution rj^c will be presented in a forthcoming publication |11) . 

This paper is organised as follows. In Section [2] we define the effective Hamiltonian 
relevant to AS" = 2 transitions. We present the details of our calculation as well as the 
analytic results in Section [3j The discussion and numerical evaluation follow in Section |4j 
In the appendix we show how our results transform under a change of the operator basis. 



^Our analysis, which uses different inputs for the physical parameters and a different error estimate, 
yields a NLO value of rjct = 0.457(73), see Sectionlij 



2 Effective Hamiltonian for Neutral Kaon Mixing 



The effective Hamiltonian Hf^J, ^ of Equation (1.5) describes the dominant contribution 



to AS* = 2 processes below the charm-quark mass scale. The loop functions 

S{x,)=x, + 0{xl), (2.1) 

_ Axt-llx1 + xl ?,xl\ogXt , . 

^^""'^ - 4(l-x,)2 " 2(1^^ ' ^^-^^ 

S{xc, xt) = -Xc log Xc + XcF{xt) + 0{xl log Xc) , (2.3) 

where the function F is defined as 



x^ — 8xt + 4 3xt 



F{xt) = * / ,, logxt - ,,, * , , (2.4) 



are used as normalisation factors of the three contributions proportional to A^, A^, and 



AcAt in Equation (1.5). In this normalisation we fix the charm-quark mass and the top- 
quark mass to rric = m^^^{mc) and rrit = mf^^iirtt) respectively in Xc and Xt- This avoids 
spurious scale dependences in rjct, Vcc, and tju, if these parameters are defined through 
Equation (|L5|f1 

2.1 The Operator Basis 

Above the charm-quark mass scale both the AS = 1 and AS = 2 effective Hamiltonians 
contribute to the Wilson coefficient of Qs2 through renormalisation group effects. In the 
following we list all operators needed for these effective Hamiltonians. They can be divided 
into three classes: Physical operators, gauge-invariant operators that vanish by the QCD 
equations of motion (EOM), and evanescent operators, that vanish algebraically in four 
space-time dimensions. 

We start with the dimension-six operators, which we choose such that problems arising 
from the 75 matrix appearing in closed fermion loops in the framework of dimensional 
regularisation do not occur |!l2|. There are two current-current operators 

Qf = {-SLl.T^qL) ® {q'LYT-d,) , 
Qf = {slI^Ql) ® {q'LYdL) , 

where g^, = |(1 — 75)5 is the left-handed chiral quark field, and q and q' are either u or c. 
The colour matrices T" are normalised such that TxT°'T'^ = S°'^/2. We use these operators 
in the linear combination 

Q'^' = ^ (1 ± ]^) Qf'^QT' = \ {{-sll.qD ® {i^l'di) ± {-sl^.qi) ® {q'irdl)) , (2.6) 
^The parameters rjcc, tju, and rjct equal rjl, 772, and rj^, respectively, as defined in Reference |7 . 
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Figure 3: Sample three-loop diagrams with IPI subdivergences that have to be subtracted by 
insertions of the EOM-vanishing operator. The corresponding IPI one- and two-loop insertions 
of Qcom are also shown. 

where a and /3 are colour indices, and Nc is the number of colours. The advantage is that 
the anomalous dimensions in the subspace of current-current operators are diagonal in this 

basiqB 

We define the QCD penguin operators as 



■^ — ^<? 
^ — ^g 
■^ — ^g 



(2.7) 



where the sum extends over the light quark fields, and we have introduced the abbreviations 

7/ilAt2/i3 ~ 7Ati7At27M3' ^^'-'• 

In order to subtract the divergences of all possible one-particle irreducible (IPI) sub- 
diagrams of the relevant Green's functions we need the following gauge-invariant EOM- 
vanishing operator 



Q 



-SLYT'^dLD^Gl. + Q^. 



eom -^Ij / -^ ^ij^ ^ av 

9 



{21 



where D^ denotes the covariant derivative, acting on the gluon field, and g^ = Anas is the 
square of the strong coupling constant. Sample diagrams are shown in Figure [3j 

The operator inducing the effective \AS\ = 2 interactions above the charm quark scale 
can be chosen as 



Qi 



m 



,.Ssi^M®{sl^^di), 



r^J' 



(2.9) 



where a and (3 again denote colour indices: Note that, according to convention, we define 
the operator with two inverse powers of the strong coupling constant in order to account 
for the logarithm already present at leading order. 



^Tliis is true beyond LO only with a suitable choice of the evanescent operators, see below. 



The use of dimensional regularisation in a theory involving fermions implies an infinite- 
dimensional Dirac algebra. In order to remove all divergences of the Green's functions, 
we have to introduce a set of evanescent operators that are non-zero in d dimensions and 
vanish algebraically in four dimensions. At the one-loop level we need 

Ef^''> = {sLl,.,,,,T^qL) ® (q'^r^^'^'T^dL) - (16 - 4e - 4e')Qf , 



Ef^'^ = {-SLl,.,,,,qL) ® {qLl'^'^'-'dL) - (16 - 4e - Ae^)Q'^'^ 



2 5 



E'i^ = {-SLl,,,,,,,,,M ® E {qi'''''''''''q) + 64g3 - 2QQ, , ^^'^^^ 

Ef^ = {-SLl,,,,,,,,,XdL) ® Y. {qi'^^'^'^'^'^T^q) + 64Q4 - 2QQ, . 
At the two-loop level, we use the following four operators: 

ET'^'^ = {sLl,,,,,,,,,,T^qL) ® {q'r^Y'^'^^''''T^dL) - ^256 - 2246 - ^6^^ Qf ' , 

Ef^'^ = (^L7..M..3.4M.?^) ® iqLl''''^''^'''''d,) - (256 - 224e - ^6^) Qf' , 

i^f = (^L7mM./^3M4M5/.6.7^L) ® V (g7^i'^^'^3/.4M5/.6Mv^) + i280Qs - 336 Q , , 

^ — ^q 

Ef^ = {-SLl,,,,,,,,,,,,,XdL) ® Y. iqY"'''"'''"''''T^q) + I28OQ4 - 336g6 • 

^ — ^q 

(2.11) 

The evanescent operators in the current-current sector are chosen such that the anoma- 
lous dimensions for the operators Q'!^ are diagonal through NNLO 13 . The remaining 



operators are chosen as in Reference 12 . 

In addition to the operator Q? with the colour structure (s27^(i2) ® (^l7'^^l)' ^^^ 
dimension-six and dimension-eight operators will also mix into an operator with the colour 
structure {s^'^^di) ® {sLl^dl). In four space-time dimensions, the latter is related to the 
former structure by a Fierz transformation. The difference of these structures is therefore 
evanescent, and correspondingly we introduce an evanescent operator of the following form: 

^^ = ^M^A) ® i4Ydl) - Q7 . (2.12) 

We choose the remaining evanescent dimension-eight operators to be 

^7'^ = -^Ml,^,.,sdl) ® {4r''''''d^,) - (16 - 4e - Ae')Q, , 
9 f^ 
2 

^s'^ = ^M^^if^^^A) ® (slr'^'^'dl) - (16 - 4e - Ae'){Q, + Ep) , 
9 r 
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,dt) ® (sf 7'^!^^'^^^*''^?) - (256 - 224e - 
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4) ® (sfT^'^^^'^^^^'^^rf^) - (256 - 224e - ^^^^^e^){Q7 + Ep) 



325 
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8 2,,2e ^ -^ '/^l'*2M3M4A»5M6M7 



4) ® (s^7'^1^2WM4/.5M6/.7^2) - 4096(^7 + Ef) . 



(2.13) 



This choice ensures that Qs2 will have the same anomalous dimension as Q+ up to NNLO. 
It is given explicitly here for the first time. 



2.2 Effective Hamiltonian 

We obtain the effective Hamiltonian valid between the electroweak and the bottom-quark 
scale by removing the top quark and the W boson as dynamical degrees of freedom from 
the Standard Model. It reads in terms of the renormalised Wilson coefficients 



njcS 
^/=5 
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lk=+- 1= 



6 



CkCiZkij + C-jZ-j-j 



Q7 + h.c. . 

(2.14) 



Here the first line represents the jAS*! = 1 part of the effective Hamiltonian, whereas 
the second line contains the jAS*! = 2 contributions. The first term in the second line is 
related to a single insertion of (5s2j induced by the top quark contribution to the Standard 
Model amplitude. The remaining terms arise from the mixing of insertions of two | AS*! = 
1 operators into the operator Q-j. The GIM mechanism leads to the absence of a A^ 
contribution to the Wilson coefficient of Q-j. The renormalisation constants Z are defined 
such that any renormalised effective amplitude, of the form Acsi = Ci{ii)Zij{ZQj)R + 
{CkCk'Zkk',i + CkZki){ZQi)R, is finite and imphcitly includes the contribution of evanescent 
operators. Here angle-brackets denote matrix elements between initial and final states i 
and /, respectively, i.e. {Qj) = {f\Qj\i). Z denotes the wave function renormalisation of 
the fields in the operator, so that {ZQi)fi are the renormalised matrix elements of the bare 
operator Q^'^'^^, where masses and gauge couplings are renormalised in the usual way. 

The effective Hamiltonian 'H'i^4 valid between the bottom- and the charm-quark scale 
looks exactly the same as "H'i^g. The only difference is induced by the presence of penguin 
operators, which explicitly depend on all light quark fields. 



Below the charm-quark scale, the charm quark is removed as a dynamical degree of 
freedom. As a consequence, the \AS\ = 1 operators can now be dropped from the effec- 
tive Lagrangian, because the matrix elements of double insertions of these operators are 
suppressed by factors of m'^/M^. The effective Hamiltonian is thus given by 



.T,|AS'|=2 _ Gj^ 



KC^^^) + KCUf^) + XcXtC^'2if^) Zs2Qs2 (2.15) 



and now only contains the \^S\ = 2 operator Qs2 defined in Equation (1.6). 



3 Calculation of rjct 

In this section we present the details of the calculation of rjct in the NNLO approximation. 
We start with the determination of the initial conditions for the Wilson coefficients at the 
electroweak scale. Afterwards we use the renormalisation group equations to evolve them 
down to the charm-quark scale, including the threshold corrections at the bottom-quark 
scale. Finally we determine the charm-top contribution to C5.2 by a matching calculation 
at the charm-quark scale. 

3.1 Initial Conditions at the Electroweak Scale 

The initial conditions for the Wilson coefficients of the dimension-six operators are available 
in the literature. In our basiSjWhere we can use a naive anticommuting 75, the results up 
to second order in ai reaq"^ 



C^{f^) = l±l[lT (11 + 6 L^) ^^ + (^ (7 ± 51) ^' T ^1 T T{x,) 
(135677 T 124095) - — (11 ^ 249) Lw + - {7 ± 51) L^ ' ' "' ^^' 
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"^Here and in the following, by the superscript in brackets we display explicitly the number of light 
quark flavours for which as is defined. 
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Figure 4: Sample two-loop Feynman diagrams contributing to the matching at the electroweak 
scale. 



CM 
CM 



An 

477 



(^^o(^.) 



10 ^^ *^ 243 






85 
162 



2 14^ 2 . 

. 81 '^ 27 ^ 



2 r r 2 

108 108 36 ^ 
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We have taken the initial conditions for C+ from Reference fl3l. The initial conditions 
for Cs . . . Ck can be found in Reference 



14 



where also the loop functions T{xt), G\{xt), 
EQ^Xt) and E\{xt) are defined. Note that in our renormalisation scheme we had to include 
an additional finite contribution for C4, as described in the appendix. We have introduced 
the abbreviation L^ = log(/i^/M^). 

With these ingredients, we can now calculate the initial conditions for the Wilson 
coefficients of the dimension-eight operators. In order to match the Green's functions in 
the Standard Model and the effective five-flavour theory, we have to compute the flnite 
parts of Feynman diagrams of the type shown in Figures [T] and |4j To this end, we perform 
a Taylor expansion in the charm-quark mass of all propagators corresponding to a charm- 
quark fleld. The constant terms cancel because of the GIM mechanism, whereas the terms 
proportional to m'^ give the leading non-vanishing contribution we are interested in. This 
procedure leads to massless vacuum integrals in the effective theory, such that only terms 
proportional to tree-level matrix elements remain. Some of these terms multiply divergent 
renormalisation constants and correspond to infrared divergences in the effective theory. 
They exactly cancel the corresponding infrared divergent terms in the Standard Model, 
leaving us with a finite result. 



Expanding the dimension-eight Wilson coefficient as 
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we obtain the following result: 

1 
2 



Cf^ (/i) = , C}') (^) = F(xt) + - - i^w 



~(2) , , 5xf - 21x,2 + 60xt - 20 , , , , 
Q H/i) = +— ^^J^ _ ^)3 ^og{xt)Lw 

12x,^-34xf-9x,3-33x,2-116xi + 36, o, , 
+ ^ 12(;-1)3 '°8 (-') 

-12xf + 27xf + 23xf + 150x2 - 108a;f + 16 , , , 
+ ^ ^. .^3 ^ ^ logout (3.3) 

6{xt - lyxt ^ ' 

-7800,T^ - 126499x3 + 191248a;2 - 129749x4 + 10400 
^ 3900(xt - l)2xi 

6x,6-llxf-8xf-29xf + 23x2-16xt + 8, , 

+ - — 3(..-i;%? ^'^<^ - ^') 



6x^ + x^ - 59x2 - 8 47x2 - 31xi + 56 ^ ^^^ 

lk?t ^' 6(xt - If 






The first line in Equation (3.3) agrees with the result obtained already by Herrlich and 
Nierste in (^ after the corresponding change of the renornialisation scheme. The two-loop 
result is new. 

3.2 Structure of the Renornialisation Group Equations 

After the determination of the initial conditions for the Wilson coefficients, the next step 
is the renornialisation group evolution to lower scales. The renornialisation group equation 
relevant for the Wilson coefficient Cy is given by: 

d ~ ~ ^ 

/i — C7(/i) = C7(/i)777 + ^ ^ Ckifi)CM%n,7 , (3.4) 

^ k=+-n=+-,3 

where 777 denotes the anomalous dimension matrix of the operator Qj, and ^kn,7 is the 
anomalous dimension tensor, describing the mixing of the dimension-six operators into Qi. 
The matrix 777 is decomposed as 777 = 752 + 27^ — 2/3, where the anomalous dimension 
of the quark mass 7^ and the /3 function are related to the factor ml/g'^ in the defini- 
tion of the operator Qj. The anomalous dimension matrix 752 is defined in terms of the 
renormalisation constants Zs2 as 

dfi 

The explicit expressions for the anomalous dimension matrix in terms of the renormalisa- 
tion constants Zs2 are given up to NNLO by 

^,(0) _ 9 7(1.1) ~ (1) _ . 7(2,1) _ ^ 7(1,1) 7(1,0) 

3.6 

r/2)_ ^7(3,1) ,7(2,1)7(1,0) ^7(1,1)7(2,0) 

752 — 0^52 " ^^S2 ^52 ~ ^^52 ^52 ' 

11 



752 = Zs2lJ^^r.^S2 ■ (3-5) 



where we only kept the non-vanishing physical contributions. Here the superscript (n, m) 
denotes the 1/e'^-pole part of the n-loop contribution. The anomalous dimension tensor is 



defined as 15 



%n,l - -{lUU'S^W + lnn'5,,,)Z,,^,,Z^l - (/i^^.n,.] Z^^ . (3.7) 

The non-vanishing contributions to the physical part of the anomalous dimension tensor 
are given in terms of the renormalisation constants by 



l'kn,l - ^^kn,l ~ '^^kn,l'^l'l' ~ 2 |Z^^', 5nn' + ^kk'Z^^, | ^fc,„,; , 

^/2) _ ^ 7(3,1) _ . A(2,l) 7(1,0) _ ^ A(l,l) 7(2,0) 
lkn,l ~ ^^kn,l ^^kn,l'^l'l ^^kn,l'^l'l 

-'^\^kk' Snn' +Skk'Z^;, j Z^.,;^,i-4^Z^.,:, 6nn' +Skk'Z^;, ^ Z^,'^, 



(3. 



(1,0) 
k'n'l 5 



where the indices fc, n and / correspond to physical operators only. 

In order to determine the renormalisation constants, we have to compute the divergent 
parts of Feynman diagrams with up to three loops, see Figure [5j We use the method 
suggested in fl6| by Chetyrkin, Misiak and Miinz for extracting the UV divergences of a 
given Feynman diagram. The renormalisation constants are then determined recursively 
by subtracting subdivergences according to Zimmermann's forest formula. As usual, we 
perform a finite renormalisation in order to ensure the vanishing of matrix elements of 
evanescent operators. An additional subtlety arises because of the presence of EOM- 
vanishing operators at second order in the effective interactions: As explained in detail in 
Reference fTT] , we have to expect non-trivial contact terms resulting from double insertions 
of Qeom and physical operators. We computed these terms explicitly, showing that non-zero 
contributions indeed occur, and subtracted them by an additional finite counterterm: 

^(2,0) ^ ^(2,0) =-( N - 1 - ^ ^-\ r3 91 

Q+Qoom,Q7 Q-Qcom,Q7 8 V '^ N^ N^ J ' 

This renormalisation ensures the validity of the equations of motion also at second order 
in the effective interactions. 

Let us now look at the equations (3.4) in more detail. It turns out that these equations 
are equivalent to the following system of eight equations [7] 

fi^D = 'y^D, (3.10) 

where the anomalous dimension matrix and the Wilson coefficients are now given by 

\ / C(^) \ 

7^ = I <7 777-7+ , D{f,) = C'+(/i)/C+(/i) , (3.11) 
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Figure 5: Sample one-, two-, and three- loop diagrams contributing to the NNLO mixing of 
dimension-six into dimension-eight operators. 



if we decompose the Wilson coefficient Cj as 



(3.12) 



This decomposition is completely arbitrary and preserved by the renormalisation group 
evolution. For instance, we may choose C^{fi) = Ci{ii) and Cf (/i) = 0. The advantage 



of (3.10) is that it has the form of a renormalisation group equation for a single operator 



insertion, and we can use the well known explicit solution (see, for instance. Reference [2]). 
We obtain the anomalous dimension matrix 7q of the operators Q_|_, Q_, Q3, . . . , Qg 
from Reference [4] by the basis transformation described in the appendix and find 
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Here and in the following, / is the number of active quark flavours. 

We denote the anomalous dimension for the double insertion of either (5+ or Q- and 
one of the operators Q+, Q--, Q3, . . . , Qe by 



7±,7 



(7±+,7, 7±-,7, 7±3,7, 7±4,7, 7±5,7, 7±6,7) 



(3.16) 



and find 
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Figure 6: Sample one-, two-, and three- loop diagrams, whose divergent parts contribute to the 
anomalous dimensions of the operator Q'j. Curly lines denote gluons, dotted lines denote ghosts, 
and solid lines denote quarks. 
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at LO, NLO, and NNLO, respectively. The LO and NLO results agree with the literature [T] 
after the corresponding change of the operator basis, described in the appendix. The NNLO 
result is new. 

In the calculation of 752 (cf. the diagrams in Figure [6]) we have chosen the evanescent 
operators in the dimension-eight sector in such a way that the anomalous dimension of 
the operator Qs2 equals the anomalous dimension of Q^ through NNLO. Consequently 
7S2 = 7+, and (l3 



,(0) 



i4" 



(2) 275 267 



9 
52 891 

675 



f-7: 



f 



260 



f 



160 



/ + 672 C: 



(3.20) 



The explicit expressions for the QCD /3 function and the anomalous dimension of the 

^^ ^ 2857 _ 5033^^ 325^, _ (3.21) 



quark mass are given by 18[j21 
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Figure 7: Feynman diagrams relevant for the threshold corrections at the bottom quark scale. 
The one-loop diagram of Qi and Q2 is the same in both theories, whereas at the two-loop level 
they receive non-trivial corrections from virtual bottom quarks. The same applies to insertions 
of the operator Qj. Because the penguin operators mix into Qeom-, we also had to calculate 
insertions of (53,...,6 with one external gluon, expanding up to the second power in the external 
momenta. 



and 

(n\ (^^ 404 40, .^^ /4432 320 \ , 280 ,, 

7i°^ = 8, ^^J,^ = — --f, j(^) = 2A98-l^-^ + —C,jf-—f. (3.22) 

3.3 Threshold Corrections at the Bottom-Quark Scale 

When we pass the bottom-quark threshold, we must perform a proper matching between 
the effective theories with five and four flavours. This threshold correction is computed by 
requiring the equality of the Green's functions in the two theories at the matching scale, 
in this case fib = 0{mb), where nib is the bottom-quark mass. 

At NNLO, there are several sources of matching corrections. The penguin operators are 
affected already at NLO, because they explicitly depend on the number of light-quark fields. 
At NNLO also the matching of the current-current and the dimension-eight operators is 
non-trivial. The source of such contributions are virtual bottom quarks in two-loop matrix 
elements of the form shown in Figure [7} In addition, also the strong coupling constant and 
the charm-quark mass are discontinuous beyond LO. 

Let us write the equality of a general amplitude in the two theories at the matching 
scale fif as 

C/-i(^/)(g/-i)(/i/) = cK^/)(g/)(/i/) , (3.23) 

the variables with subscripts / and / — 1 belonging to the /- and / — 1-fiavour theory. 
At the bottom-quark scale, we have / = 5. We parameterise the matrix elements of the 
operators as an expansion in the coupling constant defined in the corresponding /-flavour 
theory: 



{Qf)M = I 1 + ^^rf>M + (^] rfM] (0,)(«>. (3.24) 
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An additional subtlety arises, because the strong coupling constant also gets a non-trivial 
matching correction at a flavour threshold. Up to the NNLO approximation we have the 



relation 22 24 
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/ f_-t ) 
which we use to express all quantities in terms of the coupling constant aj (/W/) in the 

effective theory with / — 1 flavours. Here ruf = mf{^f) is the MS mass of the quark which 

is integrated out. Note that the matching for C7 starts at order I /as, so that by inverting 

Equation (3.25) we get a contribution already at NLO. Similarly, we need the decoupling 



relation for the charm quark mass up to NNLO 25 
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In order to display the threshold corrections explicitly, we now introduce the disconti- 
nuities 

6C^'\^,f) = CfM - Cf\M , 5rW(/x;) = rfM - rf\M , (3.27) 

of the Wilson coefficients and the matrix elements, respectively, and find for the general 
solution of Equation (3.23), in case of the dimension-six Wilson coefficients: 

5C(°)(/.^) =0, 6C^'\f,f) = -Cf\^,f)6r^'\^,f), 



SC^'\^^f) = -Cf\^^J)\5r^'\^^f) + -\og 
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Notice that the different single contributions in the last bracket may not be finite because 
of spurious IR divergences, which nevertheless cancel in the sum. The matching corrections 
look different for the dimension-eight Wilson coefficients, because of the factor l/g"^ in front 
of the operator: 
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In addition, we have to take into account the terms related to the decouphng of the charm- 
quark mass. 

At NLO, only the matrix elements of the penguin operators get non- vanishing contri- 
butions. They can be obtained from 
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where ^tq denotes the difference of the matrix elements in the subspace of dimension-six 
operators. At NNLO, we obtain the following contributions for the penguin operators: 
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where we can extract a, and a_ from 13 to find 
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(Lb = log i/J-l/fnl) here and in the following two equations). We have determined the other 
entries by calculating two-loop matrix elements of the operators Q+,Q-,Q3...6 between 
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appropriate external states (see Figure [?]), and find 
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(2)__826_128 _ 40 2 (2)__467_266 _ 23 2 
For the dimension-eight operator we find the only non- vanishing contribution 



3.4 Matching at the Charm-Quark Scale 

At the scale /Xc = 0{mc) the charm quark is removed from the theory as a dynamical degree 
of freedom, and the effective Lagrangian is now given by Equation (2.15). Requiring the 
equality of the Green's functions in both theories at the charm-quark scale leads to the 
matching condition 

6 

XI XI ^iil^c)Cj{lJ,c)ZikZji{QkQl){Hc)+C7il^c)Z77{Q7)i^c) 

i,k=+,— j,l=+,—,3 

= :^A'2if^'^)^S2{Qs2)i^^c), (3.35) 

which we use to determine the Wilson coefficient C|.2(/u)- To proceed, we parameterise the 
matrix elements in the following way: 

m = r^m^'^ , {Qs2) = rs2{Qs2)^'^ , and {Q.Q,){fi,) = '^^^r,,,s2{Qs2)^'^ ■ (3.36) 

3/7r^ 

If we take into account the explicit factor of rnl/g"^ in the definition of Qj and expand the 
Wilson coefficient Cg2 as 

^^U^) = -^.cTif^) + cTi^^) + "^CsTif^) , (3.37) 

a^if-i) 47r 

19 




Figure 8 
scale. 



Q1...6 (-/!.. .6 

: Sample one- and two-loop diagrams contributing to the matching at the charm-quark 



we find the following contributions to the matching: 
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where 6r\ (Hc) is given by Equation (3.34), with nt, and rub replaced by /Xc and rric, respec- 



tively. Notice the additional logarithms which we get by expressing al through as . 
These terms, which are numerically tiny at NLO, have been neglected in Reference [T]. 

Furthermore, we expand the charm-quark mass defined at the scale fie, viz. mc{fic)i 
about mc{mc) (see fl3]): 
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In order to evaluate the equations (3.38), (3.39), and (3.40), we have to compute the 



finite parts of one- and two-loop Feynman diagrams of the type shown in Figure pi In this 
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way we find for r±j^s2 at one loop: 
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This result agrees with the one obtained in I7j after the appropriate basis transformation. 
A two-loop matching calculation yields 
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This result is new and completes the matching onto the three-flavour theory. Now only 
a single operator contributes, and the renormalisation group evolution below the charm- 
quark scale is the same for the top-, the charm-, and the charm-top-quark contribution. 



3.5 Renormalisation Group Equations below the Charm-Quark 
Threshold 

The effective Hamiltonian valid below the charm-quark threshold contains only the single 
operator Qs2- The renormalisation group evolution is now the same for the three Wilson 
coefficients C5.2, where j = c, t, ct, and is described by the evolution matrix corresponding 
to the anomalous dimension of Qs2- 



Cs2if^) = f/(/i,/ic)<^52(/^c 



(3.45) 
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By comparing (1.5) and (2.15), we see that we can express the coefficients rjcc, f]tt, Vet as 
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The remaining /i-dependence present in (|3.46|), corresponding to the lower end of the 

(3.47) 



evolution in Equation (3.45), is absorbed into 6(/u), which equals 

6(/i) = [af (/x)]"'^"ir+(/x), 
where 
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and the exponent a+ is the so-called magic number for the operator Q^ (the magic numbers 
as well as the matrix J are defined for instance in [^). This scale dependence is cancelled 
by the corresponding scale dependence of the hadronic matrix element. 

3.6 Analytical Checks of our Calculation 

Because the calculation of the NNLO contributions to rjct is quite complex, we checked our 
results in several ways. 

First of all the calculation of the (9(100 000) Feynman diagrams as well as the renor- 
malisation, the computation of the anomalous dimensions and the matching, has been 
performed independently by the two of us, using a completely different setup of computer 
programs. On the one hand we use qgraf [26 1 for generating the diagrams; the evaluation of 



the integrals is then performed using the program packages q2e/exp/MATAD [27 28 , where 



MATAD is written in FORM 29 and based on the Integration-By-Parts algorithm 30 31 



In addition, we have written our own FORM routine in order to evaluate two-loop diagrams 
with an arbitrary number of (possibly vanishing) masses, using the algorithm described 
On the other hand, all calculations have been performed using an completely 



m 



14,32 



independent setup, based on Feynarts 33 and Mathematica 



As a check of our calculation, we verified that all anomalous dimensions, Wilson coef- 
ficients, and matrix elements are independent of the gauge-fixing parameter ^. Because of 
the complexity of the analytical expressions, for the three-loop penguin insertions we kept 
only the ffist power in ^ for our check. 

Another very useful check is the locality of the counterterms, which is an implication 
of renormalisability. In a mass independent renormalisation scheme this means that the 
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Parameter 


Value 


Ref. 


Parameter 


Value 


Ref. 


Mw 


80.399(23) GeV 
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as{Mz) 


0.1184(7) 
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163.7(1.1) GeV 
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2; 
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2] 
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17.77(12)/ps 
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^s 


1.243(28) 


9] 


Vtt 


0.5765(65) 


8' 


Bk 


0.725(26) 


9] 


Vcc 


1.43(23) 


7 


sin 2^ 


0.671(23) 


2] 









Table 1: Input parameters used in our numerical analysis. 



renormalisation factors Z depend on /i only through the coupling constants. We have 
checked this explicitly and found /i-independence of all our renormalisation constants. 

We have also checked analytically that rjct is independent of the matching scales fiw, 
fib, and fie to the considered order of the strong coupling constant, by expanding the full 
solution of the renormalisation group equations about the respective matching scale. 

As a cross-check, we confirm the NLO results of Herrlich and Nierste JTJ for the first 
time. 



4 Discussion and Numerics 

In this section we present the numerical value of r/ct at NNLO and discuss the theoretical 
uncertainty, as well as the impact on ex- Our input parameters are collected in Table [1} 

The theoretical uncertainty of r]ct is related to the truncation of the perturbation series. 
We estimate it by considering the remaining scale dependence, the different methods to 
evaluate the running strong coupling constant, and the size of the NNLO corrections. 
Varying fic from 1 to 2 GeV and fiw from 40 to 200 GeV, we find the following numerical 
value at NNLO, 



Vet 



0.496 ± 0.045^, ± 0.013^^ ± 0.002,, 



±0.001™ ±0.0002 



nit ? 



(4.1) 



where we also display the parametric uncertainties stemming from the experimental error 
on as, rric, and m^. The dependence on the scale fib is completely negligible. 

The dependence on the electroweak matching scale fiw is shown in Figure [9j We have 
plotted rjct as a function of fiw in the range from 40 GeV to 200 GeV, where we fixed the 
other scales as fib = 5 GeV and fie = 1.5 GeV, respectively. The relatively weak residual 
dependence on fi^ at NLO is slightly increased at NNLO. By contrast, the dependence 



23 



0.6 
0.5 
0.4 
0.3 
0.2 





^ 


0.5 


: 


-'' 


0.4 
0.3 
n 9 



40 80 120 160 

Mw [GeV] 



200 



Figure 9: Tjct as a function of nw at LO 
(dashed line), NLO (dashed-dotted line), and 
NNLO QCD (solid line). 
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on fifj, which is shown in Figure [TOl fixing fiw = 80 GeV and fic = 1.5 GeV, is completely 
negligible. The dependence on the scale fie is shown in Figure 11, where we vary fic 
in the range from 1 to 2 GeV, fixing fiw = 80 GeV and fib = 5 GeV. In addition, we 
have plotted rfd corresponding to three different possibilities of calculating as(/Xc) from 
the experimental input value of as{Mz)'- One method (1) is to solve the renormalisation 
group equation for as numerically. Furthermore, it is possible to compute «« by first 
determining the scale parameter Aqcd- This can be achieved by using the explicit solution 
for Aqcd without expansion in «« (method 2) or by iteratively solving this equation for 
Aqcd and from this value determining a^ (method 3). The dashed, dotted, and dashed- 



dotted lines in Figure 11, each of them representing the NLO result for rf^, correspond 
to these three possibilities of determining as, respectively. We used the mathematica 



package RunDec 37 for the numerical evaluation. Note that the difference between these 



three methods vanishes almost entirely at NNLO. On the other hand, at NLO the effect 
is sizeable and thus contributes to the theoretical uncertainty. Varying fic and fiw iii the 
same range as above, we find at NLO 



r^NLO ^ Q_457 ^ o.072„„ ± 0.01,„„ ± 0.0001„. ± 0.002^.„ ± 0.0003 



'/Jc 



f^W 



mt 1 



(4.2) 



where there error indicated by the subscript "/Xc" includes the effect of the three ways 
of determining a^. For the variation of the scale fiw we have used only method 1 for 
evaluating ag in order to avoid double-counting of the related uncertaintjQ Again we have 
included the parametric uncertainties related to a^, mc, and mt- 

The authors of Reference I7j have varied fic in the smaller range from 1.1 to 1.6 GeV, 
using a procedure equivalent to method 3 above for determining as- By looking at the 
explicit values of rjct in Figure [TT] we see that the two error bands at NLO and NNLO, 
resulting from this smaller range of fi^-, have almost no overlap. Now, with the NNLO 
results at hand, we see that our range for fi^ leads to a better estimate of the theoretical 
uncertainty. 



^Otherwise the error would amount to ±0.018 



Mw • 
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Figure 11: rjct Sts a function of fic- The LO result is represented by the double-dotted line. The 
dashed, dotted, and dashed-dotted lines correspond the NLO value of rjct, with Og evaluated by 
method 1, 2, and 3, explained in the text. The solid lines show the corresponding NNLO results; 
the ambiguity is almost cancelled. 



Looking at Figure 11 , it is striking that the scale dependence of the NLO result is barely 
reduced at NNLO. In order to understand this behaviour, let us look at the remaining 
fic dependence, which is the most pronounced, in more detail. It originates from terms 
proportional to higher powers of ag times logarithms of the renormalisation scale that are 
contained in the explicit solutions of the renormalisation group equations. These terms 
are only partially cancelled due to our truncating the perturbative expansion of the matrix 
elements at the charm-quark scale. 

We have separated the contributions to rjct of the different Wilson coefficients multiply- 



ing the matrix elements at the charm-quark scale (cf. Equations (3.38)-(3.40)). To this end 
we have chosen the operator basis as in Reference [7J, where we use the diagonal operator 
basis only in one dimension-six subspace, and Qi, • • • , Qe in the other. It turns out that 
only one contribution, proportional to the combination C_C2, shows a strong scale depen- 
dence. Although the size of the individual contributions certainly depends on the chosen 
renormalisation scheme, the general pattern is independent of this convention. It is related 
to the vanishing of the entry in the LO anomalous dimension tensor corresponding to the 
two operators Q_ and Q2- This incidence leads to a behaviour of the scale dependence for 
this single combination which would be expected from a NLO calculation, and dominates 
the scale dependence of the NNLO result. 

In general, the perturbation series for the AS" = 2 four-quark operator in an effective 
three-flavour theory is not expected to converge as well as, for instance, the perturbation 
series in Ref. 13 , where the mixing into a semileptonic operator was calculated. 
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Finally we remark that the absolute value of the NNLO correction is of the same 
order of magnitude as the range of rjct at NNLO in the interval /i^ = 1 . . . 2 GeV, so that 
using the size of the NNLO corrections as an estimate of the theoretical uncertainty yields 
approximately the same error as using the scale variation in the quoted interval. 

As a summary of the discussion above, we give the following final estimate for the 
charm-top-quark contribution to ex at NNLO: 

r/ct = 0.496 ±0.047. (4.3) 

(For comparison, an error estimate using a range for fi^ as in Reference (TI would yield 
T]ct = 0.504 ± 0.025.) The parametric uncertainty is essentially negligible with respect to 
the theoretical uncertainty. Compared to our NLO value, 

T^NLO = 0.457 ± 0.073 , (4.4) 

this corresponds to a positive shift of approximately 7%. 

Before we conclude this section, we study the impact of our calculation on the prediction 
of I ex I- To this end we use the following formulajj |38|[39] : 

|ei^| = hieCeBK\Vcb\'^\'^f]{\Vcb\^{l - p)r]uS{xt) +r]ctS{xc,Xt) - r]ccS{xc)) , (4.5) 

where 

We write f] = Rt sin/3 and 1 — p = Rt cos/3, where Rt is given by 



'Mp AMd 



Mr. V AM. 



(4.7) 



and ^s = {Fbs\/ B.s) / {Fb^v Bd) is a ratio of B meson decay constants and bag factors that 
can be computed on the lattice with high precision jol. Using the numerical values given 
in Table [l} we obtain 

\eK\ = (1.90 ± 0.04^^^ ± 0.02^,, ± 0.07^,,, ± O.IIld ± 0.22parametric) X 10-3 . (4.8) 

The first three errors correspond to rjcc, Vtt, Va, respectively. The error indicated by LD 
originates from the long-distance contribution, namely ^s, B^, and n^, which account 
for 40%, 37%, and 22% of the long-distance error, respectively. The main share of the 
parametric error stems from \Vcb\ (59%) and sin(2/3) (19%), while all other contributions 
are well below 10%. All errors have been added in quadrature. 
Compared to the prediction using the NLO value rj^^'^, 

|eNLO| ^ ^^83 ^ Q_Q4^^^ ^ 002^,, ± 0.11,,, ± O.IOld ± 0.22p,,,^etric) x 10"=^ , (4.9) 
this corresponds to a shift by approximately 3%. 



^A term proportional to RcAt/ReAc = 0(A"') has been neglected in Equation (4.5) (see Reference 38 ). 
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5 Conclusion 

We have performed a complete NNLO QCD analysis of the charm-top-quark contribution 
r]ct to the [AS"! = 2 effective Hamiltonian 'Hr^~ . We confirm the analytical results for 
r]ct obtained at NLO in Reference (TI for the first time. 

Some of our results are useful beyond rjct- The anomalous dimension of the operator Qs2 
can be employed to compute the large NNLO logarithms of 5° - -B° mixing and comprise 
part of a NNLO calculation of r]u. The NNLO matching corrections at the bottom-quark 
threshold have further applications in Kaon physics. 

Our numerical results for rjct can be summarised by a 7% positive shift in the NNLO 
prediction with respect to the NLO value, leading to rjct = 0.496 ±0.047. This corresponds 
to an enhancement of ex by roughly 3%, yielding |ex| = (1.90 ± 0.26) x 10"'^. With our 
calculation we solidified the theory prediction of ex, strengthening its role as an important 
constraint for models of new physics. 
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A Change of the Operator Basis 

In this appendix we examine how the Wilson coefficients and the anomalous dimensions 
transform under a change of the operator basis. This is important for two reasons: In 
order to find a compact form for the renormalisation group equations for double operator 
insertions, we have seen it to be useful to work in a diagonal operator basis in the subspace 
of current-current operators. However, the calculation of the dimension-six anomalous 
dimensions and Wilson coefficients has been performed in the literature in the basis given 
in [4|[l2]. Moreover, we had to transform our results in order to compare them with results 
that are available in the literature and have been calculated using yet another operator 
basis ff]. 

As is well known, a general change of the operator basis consists of a linear transfor- 
mation and a corresponding change of the renormalisation scheme ji]. Therefore, let us 
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first as a preparation derive the transformation properties of the anomalous dimensions 
for an arbitrary change of the renormahsation scheme. This generahses the aheady known 
results. Suppose we perform the following change of scheme for the Wilson coefficients 



Ci -^ a 



CjPji } 



Ck — > Ci. 



(A.l) 

'k - CjPj,, - ClCmplm,k ■ (A. 2) 

As before, we have denoted Wilson coefficients belonging to dimension-eight operators with 
a tilde and those belonging to dimension-six operators without superscript. Furthermore, 
we introduced the parameters p, p and p, which parameterise the finite transformations: 



Pij 
Pij 

Plm,k 






Pij - 






'ij 



tts .(1) 
:Plm,k + 



A-n 
A-K 

2 



(2) _ (1) (1) 
Pij Pik Pkj 



^S' 



Pik Pkj 



+ 0{al) , 



(2) 



P\'k + 0{al 



Then, in order for the effective Hamiltonian of the form 



-neff — Ci^ijQj + yCiZik + CiCjZij^kjQk 

to stay invariant, the renormahsation constants must transform as 

PikZkj 5 
Pik^kj 5 
PilPjm^lm,k ~r PilPjmPlm,pPpq^qk • 



Zij -^ Kj 

^ij ~^ ^ij 

Zij^k — ^ ^ij,A; 



(A.3) 
(A.4) 
(A.5) 

(A.6) 

(A.7) 
(A.8) 
(A.9) 



The transformation of the anomalous dimensions can now be obtained by inserting 
the transformed renormahsation constants into the defining equation for the anomalous 



dimension matrix (3.5), and the anomalous dimension tensor (3.7), respectively. In this 



way we obtain the well-known results for the case of single insertions 14,40,41 



/(O) 

'(1) 

/(2) 



7^r 



(1) 

(2) 



[p«,7(°)] 
[p(2),7(°)] 



(1) 



2/3oP 



1 
(1)1 



p«[p(S7 



) 0/(0)1 



4/3op(^) - 2Ap(i) + 2/3op(^V^'^ . 



(A.IO) 
(A.ll) 

(A.12) 



The general transformation law for the anomalous dimension tensor for double insertion 
readsO 



iij,k 
-.1(1) 
lij,k 



(0) 
nj,k ' 



lij,k + Pijjllk + ^Pii,kPo + lijAPlk 



,(0); 



(0)\ 41) 



- 1 Til ^jm + ^ilTjrn \ Pllk " 1 Pil^jm + ^ilpd \ % 



.W, 



.(1) 



■,(0) 
lm,k ■ 



(A.13) 



(A.14) 



^Note that additional finite contributions arise if we include the factor of mj./g^ in the definition of the 
dimension-eight operators. 
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Let us now examine how the anomalous dimensions and the Wilson coefficients change 
under a basis transformation. In four space-time dimensions, a change of n dimension-six 
operators Q and m dimension-eight operators Q is simply given by a linear transformation 



^i ' ^i -^ij^j 5 ^i ' ^i ^ij^ 



ij^j , 



(A.15) 



described by matrices R G GL(n), R G GL(?Ti). Under this transformation the renormali- 
sation constants change according to 



Z,^ — RikZhiR. 



'ij 



^:. 



RikZkiR. 



y 



z' 



kn,l 



Rkk'Rnn'Zk>n',l'Riii 



(A.16) 



In general, the situation is more complicated because of the presence of evanescent 
operators. As explained in detail in Reference [4] , a change of the operator basis consists 
of a linear transformation and a finite renormalisation; the latter is needed in order to 
restore the standard MS definition of the renormalisation constants. 

We can write a general transformation among all dimension-six operators as 



Q' 

E' 



R 






M 



eU + eV 



W 
1 



(A.17) 



where the matrices R and M parameterise a linear transformation among the physical and 
evanescent operators Q and E, respectively, W parameterises the addition of multiples of 
evanescent operators to the physical operators, and U and V parameterise the addition 
of multiples of e and e^ times physical operators to the evanescent operators, respectively. 
We apply a transformation of the same form to the dimension-eight operators, where 
we denote the corresponding matrices by a tilde, as before. The finite renormalisation 
constants can now be determined by requiring that an effective amplitude of the form 
CiZijiQj) + {CiZik + CiCjZij^k){Qk) be invariant under the basis transformation and be 
renormalised according to the MS prescription. 

Let us start with the anomalous dimension matrices for the mixing of dimension-six 



into dimension-six operators. The finite renormalisation induced by the change (A.17) is 
given by |4,13 



/(i,o) 

/(2,0) 



R 



WZ 



(1,0) 
EQ 



zi^i' + wz^^i^ 



1 



;7 



(°)w^U/ 



R- 



-R 



zSi^u 



zSi^v 



'(i'l)T/o,(0) 



'QE 



V-f^ 



R- 



where 



7(2,2) 



(1,1)^(1,1) 1 (0)^(1,1) 
QE ^EE ^ r, ' ^QE 



/3o^, 



(1,1) 
QE 



(A.18) 



(A.19) 



We have set W to zero in the second line of Equation (A.18) as these terms are not needed 



in our work. We now find the transformation law for the anomalous dimension matrices 
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in a straightforward manner using Equations (A. 10) to (A. 12): 



7 



'(0) 



7 



'(1) 



7 



'(2) 



R^^^^R'' - 



v'(i'O) '(0) 
^QQ '7 

^'(2,0) '(0) 
^QQ '7 



'(1,0) 



2/3o^, 



^QQ '7 



-4/3oZ 



'(2,0) 



2Az;(^°)+2/3ofz;(^'°^ 



The Wilson coefficients change according to 



^'(1,0) '(0) 

^QO '7 



7'(i.o) 
^00 



(A.20) 



c\^,) 



asjfl) /(1,0) ^ f (^s{^) \ ^/(2,0) 



^i\T 



(R-'YCif^)- 



(A.21) 



Clearly, the transformation law of the anomalous dimension matrix describing the mix- 
ing among the dimension-eight operators themselves is given by a formula completely 



analoguous to (A.20). In order to find the transformation law of the anomalous dimension 



tensor, describing the mixing of dimension-six into dimension-eight operators, and of the 
dimension-eight Wilson coefficients, we apply the same method as above. In addition to 



the finite renormalisation constants (A. 18), we now get extra finite contributions to Z: 



ij,k 



Z'^ = R,raRj Zl;^:llWu'Ul'p - Z^^^]Ul, + WmlZ^l^J + W^lZ^^''^ 



ln,p 



ml,p 



Wuzi^ja 



(1,1); 



Wniz\li^u, 



m "^mp j^pk • 



(A.22) 



Here the indices i, j, and k correspond to physical operators only. These expressions have 
never been given explicitly in the literature before. The anomalous dimension tensor then 
transforms according to 



/(O) 
lij,k 



p p „A0) p-1 



(A.23) 



'(1) 
lij,k 



,(1) 5-1 , ^/(i'°):/W^o7/(^'°)« ^ "/(°) ^/(^'°) 



(0)x , x.../(o)\ v/(^'°) /7'(i,o)x , x..v'(i,o)\v(°) 

jm I / lm,k ' 



- { lVS,m + 5ul'l 



jm I ^ lm,k 1 ^ il 'J.im "I" OaZj ,v 



(A.24) 



as can be derived easily from Equations (A. 13) and (A. 14). A special case of these formulas 



has been derived in Reference [15]. Using the definition (A.2), we see that the dimension- 
eight Wilson coefficients transform as 



CM = a{fi)R-' 



's,, + '^Mz'^m 



4tt 



- Ci{ij)Rj^^Cj{jj)R,^ 



jk 



as(/^)^/(i,o) 



An 



^mn,k 



(A.25) 
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Transformation to the Traditional Operator Basis 

The calculation of the NLO QCD corrections to rjct in [7] has been performed in a different 
basis for the physical operators than the one chosen by us. It is given by 

Q'f = (si^X) ® (Qirdi) , 

^ — '9 
Q', = isll.dl) ® Yl i^'R^'A) ' 

^ — ^Q 

Q', = {sliA) ® E ^^'r^'^r) ■ (A-26) 

Note that we have expressed the operators in terms of left- and right-handed fermion 
fields, in contrast to the definition used in ITJ. The evanescent operators chosen in ffi are 
equivalent to the following set of operators: 

E'f'^'^ = {sll,.,.,,ql) ® {qir^'^^'d^) - (16 - 46)Q'f , 



E'f^'^ = {sl^,.,.,A) ® iqir^-'^'dl) - (16 - Ae)Q"' 



2 5 



^'3'^ = K^,.,.,,dl) ® J2 i^'h^'^'^'^L) - (16 - 4e)g'3 , 

■^ — ^q 

Ef = isll,,,,,A) ® E {qiY'^'^'ql) - (16 - 46)Q'4 , 

^ — ^q 

Ef = {sll,,,.,,dl) ® Y i^'R^'"''''A) - (4 + 46)g'5 , 

^ — ^q 

Ef = {sll,,,.,A) ® E i^'Rl''''''<iR) - (4 + 4e)Q'6 ■ (A.27) 

-^ — ^q 

It turns out that in order to transform from our operator basis to the traditional one the 
following four evanescent operators must be introduced at the one-loop level in addition 
to the evanescent operators given in Equation ( 2.10[ ) (see Reference |4]): 

E'i^ = {-SLlM ® E (^7^75g) - ttQs + 77Q5 , 
^ — ^q o 6 

E^^^ = (sLl.T^dL) ® E (^T^sT'^g) - ^Q4 + ^Qe , 

^ — ^q 6 6 

4'^ = isLi,,,,,M E (^7'^'^''^'75g) - f Qs + % , 

■^ — ^q o o 

i?^'^ = {sL7,^,.,.T^dL) ® E,(^7''''''''75T'^g) - |g4 + iQe , (A.28) 

The transformation matrices R, M, W, and U representing the basis transformation 
according to Equation (A. 17), as well as the finite renormalisation induced by this trans- 
formation, can be found in I4j. The parts of the transformation matrices relevant to us are 
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given b}|j 
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, (A.29) 



(A.30) 



whereas the matrix V vanishes. They correspond to the bases 



Q ~ \Q 1 5 y 2 5 ^3' 



,Q'e), E'={E't 



/gg'(l) j^m'W 77/(1) 



^' 



E' 



(1)^ 



3 5 



eV) 



and 



Q — [Qi 5^2 )Q3)---)Q6)) E 



{Ef^'\Ef^'\Ei'\. 



^E^ 



(A.31) 



(A.32) 



in the notation of (|A.17|). The one-loop contribution to the finite renormahsation in the 



dimension-six sector is given by 
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/ 11 
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(A.33) 



The finite renormahsation relevant for the mixing of dimension-six into dimension-eight 
operators has never been calculated before. We find 



^ QQ,Qi 








-20 


20 
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20 


20 
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20 


20 


20 


20 
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3 



(A.34) 



*An additional rotation must be performed in order to change to the "diagonal" operator basis. This 
does not affect the finite renormahsation. 



32 



Transformation to the Diagonal Operator Basis 

Here we describe the change from the operator basis, where the current-current operators 
are defined as in Reference [4, 12 , to the diagonal basis, as defined in [l3] (an d in this 
work). The transformation matrices R, M, U, and V in the notation of (|A.17) are now 



given by 4, T3[ 



Ui^ 



R 




f o\ 

I 

10 

10 

10 



/I 

-1 




\0 1/ 

(z,j)G{(l,l),(2,2)} 
(z,j)G{(5,l),(6,2)} 
otherwise; 





i, 


i = J, 


M,,= 


< 20, 


(^,j)e{(9,l),(10,2)} 




lo, 


otherwise; 




r4. 


(2,j)e{(l,l),(2,2)} 


v., = < 


3712 
25 ' 

8032 
25 ' 


(^,J) = (6,2), 




lo, 


otherwise; 



(A.35) 



and the matrix W vanishes. These matrices correspond to the following bases of oper- 
ators (the roles of the primed and unprimed set of operators is reversed with respect to 
Reference fisl): 



(A.36) 



Q' = (Q+, Q_) , E' = {Ef,Ef, Ef, Ef) 



and 



Q = (Qi,Q2), E = {E['\ E','\ E['\ E','') . 



(A.37) 



All necessary renormalisation constants can be found in Reference jll. The finite renor- 
malisation is then given by 
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Here we have implicitly corrected some typos in Ref. 13 . 
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